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Abstract 

We study the (analytic) finite-size corrections in the Ising model on the strip with fixed (H — ) 
boundary conditions. We find that subdominant finite-size corrections to scaling should be to 
the form a^/N 2 ^ 1 for the free energy f N and b^ ) /N 2k - 1 for inverse correlation length £ n 1 , with 
integer value of k. We investigate the set {<ifc,&fc } by exact evaluation and their changes upon 
varying anisotropy of coupling. We find that the amplitude ratios bi^/afe remain constant upon 
varying coupling anisotropy. Such universal behavior are correctly reproduced by the conformal 
perturbative approach. 
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I. INTRODUCTION 



Finite-size scaling and corrections for critical systems have attracted much attention in 
recent decades. Although many theoretical results are now known about the critical expo- 
nents and universal relations among the leading critical amplitudes, not much information 
is available on ratios among the amplitudes in finite-size correction terms 121. New universal 

nn 

amplitude ratios have been recently presented for the Ising model [3-6]. Consider an Ising 
ferromagnet on an N x M lattice. If Ao > Ai > A2 > A3 > ... are the eigenvalues of the 
transfer matrix (TM), in the limit M — > 00 the free energy per spin, f^, and the inverse 
longitudinal spin-spin correlation length, are 

f N = -^\nA and = ^ In (A /A n ). (1) 

Here ( is a geometric factor, which is unity for the square lattice and, in triangular and 
loneycomb geometries (also for the square lattice when the TM progresses along the diagonal 
3]), corrects for the fact that the physical length added upon each application of the TM 



differs from one lattice spacing [8[. 

At the critical point T c the asymptotic finite-size scaling behavior of the critical free 
energy (Jn) and the inverse correlation lengths of an infinitely long 2D strip of finite 

width N has the form Q, Q 

lim N 2 (f N - /«,) - 2Nf surf = A, (2) 

N— >oo 

lim iVC 1 = D n , (3) 

AT->oo 

where f^, is the bulk free energy, f sur f is the surface free energy and A and D n are the 
universal constants, but may depend on the boundary conditions. In some 2D geometries, 



the values of A and D„ are known 



111 ], to be related to the conformal anomaly number 
(c), the conformal weight of the ground state (A), and the scaling dimension of the n-th 
scaling field (x n ) of the theory 

A = ^ (2I ~ A ) ' Dn = <%n ' (4) 

for strip geometry. Here ( is anisotropy parameter. The principle of unitarity of the un- 
derlying field theory restricts through the Kac formula the possible values of c and for each 
value of c only permits a finite number of possible values of A. For the 2D Ising model, we 
have c = 1/2 and the only possible values are A = 0, 1/16, 1/2. 



For Ising universality class there are five different boundary universality classes: peri- 
odic, antiperiodic, free, fixed and mixed (mixture of the last two). Two of them: periodic 
and antiperiodic boundary universal classes are in cylinder geometry and three boundary 
universal classes: free, fixed (H — ) and mixed in strip geometry. For fixed (H — ) boundary 
conditions the spins are fixed to the opposite values on two sides of the strip. The mixed 
boundary conditions corresponds to free boundary conditions on one side of the strip, and 
fixed boundary conditions on the other. In the terminology of surface critical phenomena 
these three boundary universal classes: free, fixed (H — ) and mixed correspond to "ordinary" , 
"extraordinary" and "special" surface critical behavior, respectively. 

The highest conformal weight A, and the scaling dimension x n depends on the boundary 
conditions and for fixed (H — ) boundary conditions they given by 

1 . . 

A = -, Xx = l, x 2 = 2, x 3 = 3, x n = n. (5) 



Quite recently, Izmailian and Hu [3| studied the finite size correction terms for the free 
energy per spin and the inverse correlation lengths of critical 2D Ising models on N x 00 
lattice and ID quantum Ising chain with periodic boundary conditions. They obtained 
analytic expressions for the finite-size correction coefficients a k and in the expansions 
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fc=i 

00 An) 

^ n L^i jy2fc-i ' v > 

k=l 

and find that although the finite-size correction coefficients a*., b% (for n — 1,2) themselves 
are non-universal (except for a\ and b^), the amplitude ratios for the coefficients of these 
series b k n ' /a k are universal. 

Later that result has been extended for the quantum Ising chain for antiperiodic and 
free boundary conditions 4| and for the two-dimensional (2D) Ising models on M. x 00 
lattice with the special boundary conditions studied by Brascamp and Kunz (BK) |5]. It 
was shown that Brascamp and Kunz (BK) boundary conditions [3] belong to the mixed 
boundary condition universality class although the mixed boundary condition and the BK 
boundary condition are different on one side of the long strip. 

In the present paper we will present exact calculations for a set of universal amplitude 
ratios for the Ising model with fixed (H — ) boundary conditions. We obtain analytic equa- 



1 2 L 



FIG. 1: The lattice L 

tions for and b^ (for n = 1, 2, 3, ...) in the expansions given by Eqs. © and fl?D and find 
that amplitude ratios /ojj. are universal. We will show that such universal behavior are 
correctly reproduced by the conformal perturbative approach. 



II. THE ISING MODEL ON CYLINDER 



The Ising model on a cylinder, with width N and 



leight M and free, fixed and mixed 



,1 



in terms of a double-row transfer 



boundary condition along diagonal was formulated in 
matrix D(u). They define the lattice £ as a square lattice rotated by 45 degrees, in which 
the rows have alternately L — l and L faces. Vertically the lattice have columns of V faces, 
and the periodic boundary conditions imposed in this direction by identifying the first and 
{V + l)th rows of faces. The lattice £ consists of N (zigzagging) columns edges and M 
(zigzagging) rows edges, respectively, and N and M are given by 

N = 2L M = 2L' (8) 

The lattice £ is shown in Fig. 1. We are concerned with the fixed (H — ) type of boundary 
conditions. For fixed (H — ) boundary conditions we choose the spins at the left and right 
boundaries of the lattice £ to be +1 and —1. 
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The partition function is given by 

' ~ (9) 



Znm = £exp J a^j + K ^ a i a j I 

M \ <i,j> <i,j> J 



where we define the set of spins not fixed by the boundary conditions by {a}. For fixed 
boundary conditions, the first sum is over edges in odd columns and the second sum over 
edges in even columns of the lattice. 

Consider the critical Ising model on a cylindrical lattice of N columns and M rows as 
described previously. The finite-size partition function Znm can be written as 

Z NM (u) = Tr {D(u) M ) = J2 e- 2ME "^ (10) 

n 

where the sum is over all eigenvalues of a transfer matrix D(u), written as e~ 2En<yU \ Since we 
restrict ourselves to the critical Ising model, we have sinh(2J) smh(2K) = 1. This condition 
can be conveniently parameterized by introducing a so-called spectral parameter u, so that 
sinh(2J) = cot(2w), sinh(2K) = tan(2w), with < u < 7r/4. The anisotropy parameter ( 
related to the spectral parameter u through 

C = sin4n (11) 

For isotropic system (K = J) we have u = ir/8 and C = 1. 

For fixed (H — ) boundary conditions all eigenvalues have been determined in jj], for any 
finite value of N. The energy E n associated with the eigenvalues of a transfer matrix D{u) 
are given by 

N/2 



E n {u) = In I 2~ N ' 2 f[ cosec ^ L-^ tt^ + fi k sin(4«) 



(12) 



Using the identity 

'Jfe-1/2 



k=i v 7 

the Eq. ffl~2l) can be simplified to the following form 

N/2 



E n (u) = -Jj^ln 1 + fi k sin(4w) sin ( ^- 

k=i L V + / 



(14) 



where //& = ±1. For fixed (H — ) boundary conditions should satisfied the conditions 
Tlk=i^k = — 1- This correctly yields 2 L_1 eigenvalues for the fixed (H — ) boundary cases. 



The ground state E correspond to all = 1 except /xi = —1 in Eq. (TTJJ. The exited 
levels E n are given by Eq. ffT4"]) with \i n = — 1 and all other are equal to 1: fik = 1 
(k n). Thus for the critical free energy (Nf^ = —E ) and the inverse correlation lengths 
(^jT 1 = E r — Eq) (for r = 1, 2, n) we obtain 



-i 

3 




r 1 

Sn 



(2n + 1)tt 



2(iV+l) 
where oj{x) is given by 



to 



7T 



2(iv + r 



+ w 



7T 



2(jv + i; 



cj(x) = - In [1 + sin(4w) sinx] 



The sum in Eq. (Tl5|) can be transformed as 
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E 



fc - 1/2 
iV+ 1 



7T 



1 /7T 
W — , 

2 V2/ 2 



1 " 

E 



A; + 1/2 

"ivTT 



7T 



(2n + 1)tt 
2(iV + l) 



(19) 



(20) 



(21) 

l^ . Suppose that 



fc=l x ' k=0 

and can be handled by using the Euler-Maclaurin summation formula 
F(x) together with its derivatives is continuous within the interval (a,b). Then the general 
Euler-Maclaurin summation formula states 



N-l 



1 fo _ _ ^h*- 1 



F(a + nh + ah) = - I F(r) dr + ^ — B,(a) (F^Q) - F^a)) (22) 



fc! 



n=0 fc=l 

where < a < 1, h = (b — a)/N and Bfe(a) are so-called Bernoulli polynomials defined in 
terms of the Bernoulli numbers B p by 

B k (a)=J2 B P71—^ ak ' P ( 23 ) 



p=0 



' (fc — p)\p\ 



Indeed, B n (0) = B n . Bernoulli polynomials satisfy the identity: 

B n {l/2) = (2 1 -™ - 1) B n 



(24) 
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Using the Euler-Maclaurin summation formula given by Eq. ( 122]) with a = 0, b = ir, a = 1/2 
and F(x) = u(x) the asymptotic expansion of the critical free energy can be written in 
the following form 

= ^~t^{^^) feF< (25) 



i (2P-1)! 
23tt sin 4m 247(sin 3 4m - \ sin 4m) 



7T 



3 



^" r/ 48(iV+l) 5760 \N+l) + '"' 

where 

1 /' 7r 

foo = - uj(x)dx (26) 

1 1 r 

fsurf = -7 In (1 + sin 4m) + - / uj(x)dx (27) 
4 ^ Jo 

and Afc is the coefficients in the Taylor expansion of the u;(x): 

00 . 

p=0 y 

with Ao = 0, Ai = 2 sin 4m, A 2 = — | sin4u 2 , A 3 = sin 3 4m — | sin4u, .... 

Using the Taylor expansion of the u(x) the asymptotic expansion of the critical inverse 
correlation lengths i^T 1 for r = 1, 2, 3, n can be written in the following form 

r 1 - f 2(3 *~ 1 - 1)A2 - 1 2P " (29) 

JV + 1 + 12 ^ivTT 1 + "' 

r 1 - V 2(52P " ~ ^ 2P_1 (so 

— h r) +..-, 



-f 2 2 P~ 1 (2p-l)! V^+i 
7rsin4M 13 (sin 3 4m — | sin 4m) / jv 



^ 2 2 P- 1 (2p-l)\ \N+1 



27rsin4M 31 (sin 3 4m — \ sin 4m) / ;;r 



N + l 6 \N + 1 

00 2(7 2 ^ - l)A 2p _i ' - ^ 2p ^ 



4i 2 2 f" 1 (2p-l)! \N+lJ 



37rsin4M 57(sin 3 4m — | sin 4m) / 7r 



3 



iV + 1 + 4 liV + 1 ' 



00 2[(2n + l) 2 "- 1 -l]A 2p _ 1 / vr ^ 2p " ] 



2 2 p-!(2 P -1)! P ^iV + lJ ' ^ 

M7r sin 4m [(2n+ l) 3 - 1] (sin 3 4m- ± sin 4m) / 7r 
iV + 1 + 24 ViV + 1 1 : 
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Eqs. f[2"5"j) and f[2"P"|) - ( 13"2"j) imply that the ratios of amplitudes of the (N + 1) _ ( 2 p _1 ) correction 
terms in the free energy and the inverse correlation lengths expansions, i.e., b p /a p should 
not depend on the spectral parameter u and given by 

b p n) (2n + l) 2 f- 1 -l 



a - 2**-*B 2p (l/2) 1 

p 2p ~ x 



(33) 



For p = 1 we have 



= n 34 

ai 23 v 7 



and for p = 2 we have 



fel n) 240 r/ . o -i 480 . o 

^7 = ~247 K 2n + 1} " 1] = ~247 n(4n + 6n + 3) ' (35) 



III. PERTURBATED CONFORMAL FIELD THEORY 

The finite-size corrections to Eqs. (E D and (jSJ) can be calculated by the means of a 



perturbated conformal field theory 
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151 ] . In general, any lattice Hamiltonian will contain 



correction terms to the critical Hamiltonian H c 



,JV/2 

H = H c +y]g p (j> p (v)dv, (36) 

„ J —N/2 



P 

where p p is a non- universal constant and <p p (v) is a perturbative conformal field. Below 
we will consider the case with only one perturbative conformal field, say 4>i(v). Then the 
eigenvalues of H are 

N/2 
-N/2 

where E n>c are the critical eigenvalues of H . The matrix element < n\4>iiy)\n > can be com- 
muted in terms of the universal structure constants (C n / n ) of the operator product expansion 
14) ]: < n\(j)i(v)\n >= (2tc/N) Xl C n i n , where X\ is the scaling dimension of the conformal field 

4>i{v). The energy gaps (E n — E Q ) and the ground-state energy (E ) can be written as 

2vr ,„ „ s (2ix\ Xl ~ x 



,N/2 

E n = E niC + gi <n\(j)i{v)\n> dv + ..., (37) 

J -N/2 



E n -E = —x n + 2ir 9l (C nln - C olQ ) l-j + . . . , (38) 
E = E 0)C + 2ix 9l C m - + . . . . (39) 



Note, that the ground state energy Eq and the energy gaps (E n — Eq) of a quantum spin 
chain are, respectively, the quantum analogies of the free energy and inverse correlation 
lengths for the Ising model; that is, 

Nf N ^-E , and ^ 1 ^E n -E . (40) 



For the 2D Ising model, one finds 16( that the leading finite-size corrections (1/iV 3 ) can be 
described by the Hamiltonian given by Eq. (136]) with a single perturbative conformal field 
4>i(v) = L 2 _ 2 (v) with scaling dimension xi — 4 . 

In order to obtain the corrections we need the matrix elements < n\L 2 _ 2 {y)\n >, which 



have already been computed by Reinicke 



< A + r\L 2 _JA + r > 



N 



2vr\ 4 f 49 / 5 r(2A + r)(5A + l) 



+ (A + r) A--+ A ■ ■ ' (11) 



11520 v 7 V 24 (A + 1)(2A + 1) 
The universal structure constants C212, Cm and Cqiq can be obtained from the matrix element 
< n\L 2 _ 2 (v)\n >= (2%/N) Xl C n i n , where x\ = 4 is the scaling dimension of the conformal field 
Ll 2 (y). 

At the critical point the spectra of the Hamiltonian with fixed (H — ) boundary conditions 
can be understood in terms of irreducible representations A of a single Virasoro algebra with 
values of A is |. 

For fixed (H — ) boundary conditions the ground state |0 > and the excited states |r > 
(r = 1, 2, n) are given by 

\r > = \A = -,r = k>, for k = l,2,...,n (42) 

After reaching this point, one can easily compute the universal structure constants C r \ r 

(r = 0, 1, 2, ...,n) for fixed (H — ) boundary conditions. The values of C r \ r can be obtained 

from Eqs. (j4ip and fj4"2l and given by: 

1729 7r(4r 2 + 6r + 3) 
rlr ~ 11520 + 24 
Cqio = 1729/11520, C m = 45409/11520, C 2n = 210049/11520, .... (43) 

Equations (!38|) and (1391) implies that the ratios of first-order corrections amplitudes for 
E n — E (^n 1 ) and — E (f N ) is universal and equal to (C i O — C n [ n )/C O io, which is consistent 
with Eq. d3SD 

/o\ Cqio — C n i n 480 , 2 . , . 

rn(2) = ^ = -— n(4n + 6n + 3) (44) 
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IV. SUMMARY 



In this paper we present exact calculations for all coefficients in the asymptotic expansions 
given by Eqs. and (J7J) for an infinite number of specific levels E n (n = 0,1,2,...), for 
which we find that the ratios b p n ^ /a p are indeed universal and doesn't depend on the spectral 
parameter u. We find that such universal behavior are correctly reproduced by the conformal 
perturbative approach. 
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